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Abstract 

A theoretical framework is developped leading to a sound derivation of Periodic Boundary Conditions 
(PBCs) for the analysis of domains smaller then the Unit Cells (UCs), named reduced Unit Cells 
(rUCs), by exploiting non-orthogonal translations and symmetries. A particular type of UCs, Offset- 
reduced Unit Cells (OrUCs) are highlighted. These enable the reduction of the analysis domain of 
the traditionally defined UCs without any loading restriction. The relevance of the framework and 
its application to any periodic structure is illustrated through two practical examples: 3D woven and 
honeycomb. 

1. Introduction 

Numerical analysis of periodic materials and structures has proven to be an extremely powerful tool. 
It provides detailed information, such as failure initiation sites and stress-strain at smaller scales 
(meso/micro) It has been successfully used to determine homogenised properties, study the detailed 
stress-strain fields at nano- and microscopic scales to obtain structural damage initiation conditions 
and sites, as well as to simulate damage development and associated deterioration of the 
homogenised mechanical properties [Ij. Several works can be found discussing the application of 
periodic boundary conditions to representative regions, e.g. [2-5J. For periodic structures, the Unit 
Cell (UC) is used as the representative region, and the analysis is performed by applying periodic 
displacement boundary conditions. The topological complexity of many UCs found in practice, such 
as in typical woven composites, often leads to unpractical modelling and analysis times. For this 
reason, internal symmetries of the UCs must whenever possible be exploited to reduce the analysis 
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domain further (provided the appropriate boundary conditions are apphed), thus reducing both 
modelling and analysis time. 

A comprehensive study on the determination of reduced Unit Cells (rlJCs) for UD and particle 
reinforced composites was performed by Li [6l|7] and Li and Wongsto [8j. Different rlJCs, loading 
cases and correspondent boundary conditions were determined and presented in detail. Applied to 
textile composites, Tang and Whitcomb ^ proposed a general framework for determining rUCs. 

In the first part of the present work, the derivation of the framework proposed in [9] is revisited and 
some of its building blocks redefined, resulting in a different, formally defined and more concise 
formulation. The framework proposed by Tang and Whitcomb [9j requires the distinction of two 
different cases of equivalence between subcells: (i) equivalence is obtained by a symmetry operation 
or a translation, and (ii) equivalence is obtained by the combination of a symmetry operation and a 
translation. In the second case an additional vector of constants r (see Tang and Whitcomb [9j) 
needs to be considered when applying the boundary conditions. The non-zero components of this 
vector are tabulated for different cases and are determined by the FEM as part of the solution. The 
formulation derived in the present work is more generic, in that no cases need to be treated 
separately, and mathematically complete in that no vector r needs to be determined from tabulated 
data. All terms in the equation that assigns the periodic boundary conditions for a rlJC are fully 
defined, simplifying the formulation and consequently their use. 

In the second part of this paper, the application of the formulation developed and its potential 
is illustrated through two practical examples: 3D woven composites and honeycombs. Additionally, 
particular attention is given to Offset-reduced Unit Cells as they allow the domain reduction without 
load restrictions. 

2. Equivalence framework 

In this section, the equivalence framework is formally defined. It is based on four concepts: physical 
equivalence, load equivalence, periodicity and load admissibility. In the following sections each of 
these concepts is detailed. 
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2.1. Physical equivalence 

Consider a domain D in space and within it a sub-domain E. The latter has a defined boundary, 
Local Coordinate System (LCS) Oexi/z, and a certain spatial distribution of n physical properties P* 
with i G {1, n} . Each of these physical properties P* are expressed as a tensor written in the LCS 
ofE, i.e. P^^;. 

Definition 1. Two distinct sub- domains E and E are Physically Equivalent, denoted: 

E=E (1) 

if for every point A in E there is a point A in E such that, for each physical property i, 

= x| A P'^; {A) = V% (A) (2) 

and vice-versa. 

In Eq. jij x;^ and x-| are the coordinate vectors of A and A given in the LCS Osxyz and O^xyz 
associated with E and E, respectively, Fig. [l] The points A and A for which Eq. [l]is verified are 
designated as physically equivalent points. 




Figure 1: Geometrical relation between equivalent points. 
2.2. Periodicity and Unit Cell 

Across the literature, different definitions can be found for periodic structure and UC. In the present 
work, periodic structure and UC are defined based on the concept of physical equivalence. 
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Definition 2. A domain D is periodic if it can be reconstructed by tessellation of, non- overlapping, 
physically equivalent sub- domains Ei with parallel LCS, i.e. if for all i ^ j : 

Ei=Ej A OE^xyz \\ Oe^xyz (3) 

The smallest sub-domain verifying the periodicity definition is designated as an Unit Cell. 

2.3. Loading equivalence 

The concept of load equivalence (see Tang and Whitcomb [9]) provides a relation between 
physically-equivalent sub-domains, once the structure they are part of is loaded. Let us consider a 
periodic structure as defined in the previous section. 

Definition 3. Load equivalence between two physically equivalent points A and A is verified if the 
strains and stresses at these points, given in the LCS of the sub- domains, can be related by: 

EE {A) = (A) (4) 

aE {A) = (i) , (5) 

where the load reversal factor, 7 = ±1^ is used to enforce the equivalence between fields of physically 
equivalent sub-domains. 

For Eqs. [4] and [5] to hold, the length scale of the loading variation must be larger than the length 
scale of the sub-domains, such that an approximately periodic variation of the strains and stress fields 
is assured. If a structure is entirely composed by load equivalent sub-domains, its response can be 
obtained by analysing one of these domains alone, instead of analysing the entire structure. However, 
the appropriate boundary conditions have to be applied. These guarantee that the sub-domain, 
although isolated, has the same response it would have if it was embedded in the structure. 

2.4- Sub-domain admissibility 

Not all physically equivalent sub-domains can be used to analyse the response of a periodic structure 
under all loading conditions. The use of sub-domains smaller than the UC to analyse the response of 
a periodic structure is restricted by the relations between the LCS of these sub-domains. The 
sufficient and necessary condition for admissibility of a sub-domain to be used in the analysis of a 
periodic structure is derived below. 
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Average and fluctuation fields 

For convenience, the strain field of a sub-domain at a point A is decomposed as the sum of a volume 
average and a fluctuation term, see Figj2j 

e{A) = {s)+e*iA), (6) 

where (•) = y fy •dV is the volume average operator over the volume and e* is the fluctuation 
term, see Suquet [lOj. It is possible to find the displacements at a given point by integration of Eq. |6] 
Assuming small displacements and no rigid body rotations, the displacement relative to the origin of 
a LCS, attached to the subdomain, comes as: 

u(A) = (e)-x^ + u*(A). (7) 




X 



Figure 2: Idealised relation between the fluctuation and average fields of strain and displacement in a 
periodic structure. 

Relations between fields of two equivalent points in a global coordinate system 

Knowing that the coordinate vectors of two equivalent points A and A given in their LCSs are 

identical, Eq. [2j they can be related in the LCS of the sub-domain E by: 

x^ = T(xi-x°*) (8) 
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where T is the transformation matrix between the LCSs of E and and x^^ is the position vector 
of the origin of the LCS of the sub-domain E given in the LCS of the sub-domain E^ Fig. [l] 
Similarly, using Eq. [4j the strains at two equivalent points can be related in the LCS of E by: 

se {A) = ^Tse (i) T*. (9) 

The relation between the volume average of the strain of the equivalent sub-domains E and E^ in the 
LCS of the first, can be obtained directly by integrating Eq. [9j 

(£)| = 7T{£)|t*, (10) 

where the lower index of (•) refers to the coordinate system, and the upper index to the domain over 
which the volume average was taken. Decomposing the strain field in Eq. [9] into its average and 
fluctuation parts and using Eq. [Toj the relation between the strain fluctuations field of two equivalent 
points is obtained: 

e*E {A) = ^Tel (i) T*. (11) 
In general, the displacement of a point A can be obtained from: 

A A A A 

u(x^) = uo + / du = uo^ V.U(ix = uq + / ed:>c + / l^dx, (12) 

where £ = ^ (V.u + V.u^), ^ = ^ (V.u — V.u^) and x^ is the coordinate vector of point A. 
Considering the structure has no rigid body motion uq = , nor rotation 1] = 0, the displacement of 
a point A can be simply obtained by: 



u(x^) = / 

^x=0 



ed^. (13) 



The strain fluctuations of two equivalent points belonging to the sub-domains E and E are related in 
the LCS of E by (see also Eq. Ill]): 



(xi)=7Ts^(xi)T*. (14) 
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Knowing that two equivalent points are related by (see also Eq. [8|: 

and that the displacement at the origin is equal to zero, integrating Eq. [l4]it is possible to obtain: 

A A 

[ ^ ef (x) dx = / ^ 7Te|* (t*xb + x^^) T*Tdxs 

Jxe=0 Jxe=0 ^ ' 

ul;(xi) = 7T {u^(T*xi + x^-) - ul;(x^-)} . (15) 



Equation 15 provides a relation between the displacement perturbations of two equivalent points 



given in the LCS of one of the sub-domains. Apart from u*(x^^), all variables are known; below it is 
shown that u* [y^^^ = 0. 



According to the definition of periodicity, a periodic structure can be reconstructed from physically 
equivalent sub-domains with parallel coordinate systems. The strain fields at two equivalent points 
belonging to different sub-domains are related by: 



(X^) 



7T^^(T^x^+x^^)T^ (16) 



If we consider that the sub-domains are UCs, since the coordinate systems are parallel, the matrix T 
will be equal to the identity matrix. Moreover, all equivalent sub-domains will be admissible and have 



a load reversal factor 7 = 1. Equation [16] can then be simply written as: 

£^(xi)=£^;(x^ + X°-). (17) 

If the two UCs being considered are adjacent, a vector d = x^^ can be defined and Eq. [Tt] can be 
generalized for any point x of the structure: 

£(x) =£(x + d), (18) 

where d is commonly named the periodicity vector, and corresponds to the period of the function s (x) 
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The integral of a periodic function / of period D can always be written as: 

f{t)dt = g{t) + ft + C, (19) 
where g {t) is also a periodic function of period f is the average the periodic function /, and C is 



/■ 



a constant [llj. Using Eq. 19 it is possible to write: 

u* (x) = ^* (x) dx = h (x) . (20) 

The average term that would appear in Eq. |20] is zero since by definition (x) has zero average. 
Additionally, knowing that at the origin u* (x) is equal to zero, one can conlude that C will also be 
zero and thus u* (x) will be a periodic function with zero average. Using the above result, and knowing 
that the integration over a period of a periodic function with zero average is equal to zero, one can 
integrate both sides of Eq. [15] over a period: 

y^"^ u*(x)dx = 7 J"^"^ T {u*(T^x + x^^) - u*(x^^)} ^ix 

= 7 r^^u*(x^-)dx, (21) 

J X 

obtaining: 

u*(x^-) = 0. (22) 

Substituting Eq. [22] in Eq. [Tsj the relation between the displacement perturbations at two equivalent 
points can be finally obtained: 

uU^i)=^Tu%{4). (23) 

Evaluation of the sub- domain admissibility 

For a sub-domain to be admissible, the volume average (homogenised) strain calculated for this 
sub-domain on a given reference system must equal that volume average on any other sub-domain (on 
the same reference system), as the volume average is a homogenised entity, hence independent of the 
sub-domain where it was calculated. From load equivalence, the strains at physically-equivalent points 
are related (Eq. 5). Eq. 11 is obtained by simply integrating this relation over the sub-domain, but 
does not enforce directly that the volume average strain is a macroscopic entity independent of the 
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particular sub-domain. For the sub-domain to be admissible, the following condition must be verified: 

{efs = i^fE (24) 
as, for a sub-domain to be admissible, the homogenised strain on a given reference system (in this 
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case E) must be the same for any sub-domain (in this case E and E). Therefore, Eq. 10 with Eq. 
lead to the condition of sub-domain admissibility, as defined below. 

Definition 4. A given sub-domain E is admissible for the analysis of a periodic structure under a 
given loading {s)^, ifTi and 7^ correspondent to any other sub-domain Ei are such that, for all Ei'. 

(s)|=7.T,(£)f T* (25) 

Equation [25] can be used to, for a given loading, determine the load reversal factors 7^ associated 
with each of the sub-domains. The admissibility of a subdomain for structural analysis leads to the 
definition of a rlJC. 

Definition 5. A reduced Unit Cell is a domain, smaller than the Unit Cell, that can be used to 
determine the response of a periodic structure to a given loading. The condition to be verified by a 
reduced Unit Cell in structural analysis is defined by Eq. 

3. Derivation of Periodic Boundary conditions 

To ensure the response of a periodic structure under a given loading can be determined from the 
response of a rlJC, the appropriate boundary conditions that must applied to the latter need to be 
determined. In this section, the equivalence framework, presented previously, is used to derive the 
periodic boundary conditions for the analysis of rlJCs. 

Consider two adjacent sub-domains E and E that are physically and load equivalent. If a point A 
belonging to E is chosen to be at the boundary of the sub-domain E^ then its equivalent point A is 
also be at the boundary of E. Since both points A and A belong to E^ the displacement at each 
point can be obtained using Eq. [7| 

u{A) = (£)x^ + u*(A) (26) 
u(i) = (£)x^ + u*(i) (27) 
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All quantities in Eqs. [26] and 27 are written in the LCS of and the volume average is taken over 
the sub-domain E (the subscript will be omitted hereafter for convenience). Since both points are 
equivalent, their positions are related by Eq. [8] leading to: 

u {A) = (e) T (x^ - x^^) + u* {A) (28) 

Knowing that the displacement fluctuations at two equivalent points are related by Eq. [23j if Eq. [27| 
is multiplied by 7T and then subtracted to Eq. [26j the displacement fluctuations cancel, leading to: 

u {A) - 7TU (i) = {{e) T - 7T {e)) x^ - {e) Tx^^ (29) 

Provided the sub-domain E is admissible, see Definition 4, the term ((e) T — 7T (e)) is zero. Using 
this result, Eq. [29] can be simplified to Eq. [30j which is the main outcome of this analysis and can be 
used directly to apply periodic boundary conditions to a sub-domain: 

u(A)-7Tu(i) = -(£)Tx^^ (30) 

Once a displacement constraint equation is associated to all points at the boundary of the 
sub-domain E, loading can be applied by prescribing a volume average strain (e). 

It is relevant to notice that the displacement constraint equation traditionally used to impose 



periodic boundary conditions on a UC, see Suquet [lOj for example, is a particular case of Eq. 30 



where the matrix T is equal to the identity matrix I, since the LCSs of the UCs are parallel by 
definition and consequently, from the sub-domain admissibility evaluation, the load reversal factor is 
equal to one. It is important to highlight the differences between the result obtained above and the 



one obtained in Tang and Whitcomb [9j; as referred before, Eq. 30 is completely generic and self 
sufficient: no distinction needs to be made, in the current formulation, between the type of operation 
needed to achieve equivalence between subdomains. Moreover, all terms in Eq. [30]are fully defined, 
and can therefore be readily used to prescribe periodic boundary conditions to a given subdomain. 
3.1. Offset reduced Unit Cells 

According to the periodicity definition given in |2.2[ a UC is the smallest sub-domain that allows a 
periodic structure to be reconstructed by tessellation of sub-domains that are physically equivalent to 
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the UC and have parahel LCS. Nevertheless, in most apphcations the UC is defined such that the 



LCS are not only parallel but orthogonally translated, Fig. 3a However, smaller UCs can in general 
be defined if non-orthogonal translations are considered. Fig. [3bj Although, according to the 
definition, the representative sub-domains obtained through non-orthogonal translation are in fact 
UCs, in the present paper they are referred to as Offset-reduced Unit Cells, since they lead to a 
reduction in the domain of the traditionally defined UCs, Fig. |3] 

An important feature of OrUCs is that all loading combinations are admissible. This key feature has 
been identified by Li [7J and used in the derivation of PBCs for rUCs of UD composites, and cracked 
laminates [12j. It is relevant to highlight that using the present formulation this feature comes as a 
natural result: since the LCS of all sub-domains are parallel, they relate to each other by the identity 



matrix, i.e. T = I, as a consequence Eq. [25] is always verified and therefore all loading cases are 
admissible. 





(a) Unit Cell (b) Offset reduced Unit 

Cell 

Figure 3: a) Unit Ceh (UC) and b) Offset reduced Unit Ceh (OrUC) of a 2 x 2 Twih weave 
4. Applications 

In the present section two applications of the formulation presented previously are illustrated. 
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4.1. 3D Woven Composites 

The UCs of 3D woven composites can be significantly larger than their 2D counterparts, mostly due 
to the more intricate reinforcement architecture and 3D nature. Therefore, the domain reduction 
enabled by the use of rlJCs can be very significant. Figure |4j shows an UC, an OrlJC and a rlJC of a 
given 3D woven architecture, highlighting the domain reduction achieved: OrlJC and rlJC reduce the 
analysis domain to Y? and 1/28 of the UC, respectively. 

To define the periodic boundary conditions for the analysis of the rlJC of Fig. |4j Eq. [30j the 
geometric relations between equivalent points at the rlJC boundary need to be found. These are 
obtained by applying Eq. [8] to the equivalent domains at the boundary of the rlJC and are given in 
Table [1] and illustrated in Fig. [5] Since in general, T 7^ I, the load admissibility needs to be evaluated 
and 7, for each adjacent subdomain, determined. This is performed evaluating Eq. |25j and summarized 
in Table [2] Given a certain loading and using the data from Tables [T] and |2j Eq. |30]can be fully defined 
and periodic boundary conditions prescribed to the rUC. 
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FACES 




FACE 6 

Figure 5: Geometrical relations between equivalent points at the boundary of the rUC. 



Table 1: Geometrical relations between equivalent points at the boundary for the 3D woven rUC. /, w 
and t are respectively, the length width and and thickness of the rlJC. 
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E2 



E3 



E4 



E5 



Ee 





. = f 

-i < y < 
L 2 - ^ - 2 J 




x=f 
<V < h 
L 2 - ^ - 2 J 




y=i 

L 2-^-2 J 




. = -f 

< y < i 
L 2 - ^ - 2 J 




-i < y < 

2-^-2 J 




-f < ^ ^ f 

^ = -i 

_± < 2 < ± 

2-^-2 J 






A 

W — X-^ 

-I - x^ 
2 ,2 




A 

x-^ — w 

x^ - ^ 
2 ^ 2 

^3 




-X^ + I 

^3 




A 

— x-^ — W 
^2 J 2 
^3 




+ w 
^2+2 
^3 
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Table 2: Admissible loading cases and respective value of the load reversal factor 7^, correspondent to 
each adjacent sub-domain Ei^ for the 3D woven rlJC. 



7i Admissible loading 



CTll 


O- 12 







<722 











<733 








<7 13 








<723 


0-31 


0-32 






4.2. Honeycombs 

Honeycombs are other example of an extensively used periodic structure for which UC modelling and 
analysis can be simplified by the use of rlJCs. Figure [6] shows the UC and rlJC for a honeycomb 
structure. Following the procedure described previously, the geometrical relations between equivalent 
points at the boundary are first determined, Table [3j and Figs ?? and[7a[ The load admissibility is 
evaluated and the load reversal factors 7^ found. Table [4j 




Figure 6: Unit cell (UC) and reduced Unit Cell (rUC) of a honeycomb structure. 
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(a) rUC and adjacent sub-domains (b) Geometrical relations between equival- 

ent points at the rUC boundary 

Figure 7: 



Table 3: Geometrical relations between equivalent points at the boundary of the honeycomb rlJC. 
w are respectively, the length and width of the rlJC. 



El E2 E3 



T 


■ -1 " 
-1 




" -1 ■ 
1 




" 1 ■ 
-1 






w 





—w 







-I 






_1 < 7/^< 1 
2 — — 2 




2 — y — 2 




2 < < 2 





-x^ + w 




A 

— x-^ — w 




i 




















— X2 — I 



Table 4: Admissible loading cases and respective value of the load reversal factor 7^, correspondent to 
each adjacent sub-domain Ei^ for the rlJC of the honeycomb. 



7i Admissible loading 


Case 1 [111] 


o-ii 





<722 




Case 2 [ 1 -1 -1 ] 




0-21 


<7 12 
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5. Conclusions 



A theoretical framework leading to a sound derivation of PBCs for the analysis of domains smaller 
then the Unit Cells (UCs), named reduced Unit Cells (rUCs), by exploiting non-orthogonal 
translations and symmetries whithin the UC was developped. The investment in defining the 
problem formally resulted in a simple and readily usable formulation. The method is applied to two 
different periodic structures illustrating the potential of the rUC concept. Offset reduced Unit Cells 
are highlighted as a particular case with interesting features, allowing the analysis of domains smaller 
than the UC without any load restrictions. 
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